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Abstract 

We calculate the long-distance matrix elements for the decays of the lowest- 
lying S- and P-wave states of charmonium and bottomonium in quenched 
lattice QCD, using a nonrelativistic formulation for the heavy quarks. (The 
short-distance coefficients are known from perturbation theory.) In particular, 
we present the first calculation from QCD first principles of the color-octet 
contribution to P-wave decay — a contribution that is absent in potential mod- 
els. We also give the relations between the lattice matrix elements and their 
continuum counterparts through one-loop order in perturbation theory. 
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Typeset using REVTgX 



Heavy quarkonium systems (charmonium, bottomonium) are nonrelativistic: in the 
CM frame, the average quark velocity v satisfies w 2 < 1 (v 2 w 0.3 for charmonium, and 
v 2 ~ 0.1 for bottomonium). Bodwin, Braaten and Lepage (BBL) [TJ have shown, within 
the framework of nonrelativistic QCD (NRQCD), that the smallness of v 2 allows one to 
express a quarkonium decay rate as a sum of terms, each of which consists of a long- 
distance (distance ~ 1/Mqv) matrix element of a 4-fermion operator in the quarkonium 
state multiplied by a short- distance (distance ~ 1/Mq) parton-level decay rate, which may 
be calculated perturbatively. In particular, the decay rates for S-wave quarkonium through 
next-to-leading order in v 2 are given by 

T( 2s+1 Sj - X) = g 1 ( 2s+1 Sj)2lmf 1 C s+1 S J )/M 2 

+ F 1 (S)2Img 1 ( 2s+1 S J )/M*. (1) 

To the lowest non-trivial order in v 2 , the P-wave decay rate is given by 

TC^Pj -> X) = H l (P)2\mf 1 ( 2s+1 Pj)/M A Q 

+ Hs(P)2Imf s e s+1 Sj)/M 2 . (2) 

The /'s and g's are proportional to the short- distance rates for the annihilation of a QQ 
pair from the indicated 2s+1 Lj state, while Qi, J-\, Tlx, and Tis are the long-distance matrix 
elements. The subscripts 1 and 8 indicate whether the QQ pair is in a relative color-singlet 
or color-octet state. In this paper, we report a lattice calculation of the long-distance matrix 
elements in QCD for the lowest-lying S- and P-wave charmonium and bottomonium states. 
The calculation of H,$ yields the first result for a heavy-quark color-octet matrix element 
that is based on QCD first principles. 

The long-distance matrix elements are defined by 

Si = fSltfxxW, (3a) 
J- 1 =( 1 5|^t ( _!5)2 x |i 5); (3b) 

Hx = ( x P\1>*(i/2) D X.X f (i/2) D ^P), (3c) 
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n 8 = CPi^T-xx^yfp). 



(3d) 



The terms proportional to Q\ and Tti in the decay rates are those that appear in the 
conventional, color-singlet model ||. In the vacuum-saturation approximation which is 
correct up to terms of order t> 4 , Q\ = ^r|-R,s(0)| 2 and Hi = ^:|-Rp(0)| 2 , where -R(O) is the 
radial wavefunction at the origin and -R'(O) is the derivative of the radial wavefunction at the 
origin. The matrix-element forms of @ serve to define a regularized -R(O) and a regularized 
R'(0) in QCD. 

In contrast, the term in the P-wave decay rate that is proportional to Tig is absent 
in the color-singlet model. Tig is the probability to find a QQg component in P-wave 
quarkonium, with the QQ pair in a relative S-wave, color-octet state. As such, it corresponds 
to a true field-theoretic effect of QCD that is absent in any potential model of quarkonium. 



FIG. 1. Lattice calculation of matrix element of 4-fermion operator. The large discs represent 
the sources and sinks; the smaller discs represent the 4-fermion and point source operators. The 
lines are the nonrelativistic quark propagators. 

On the lattice, the long-distance matrix elements are obtained from the graphs of 
Fig. |l|. The upper and lower graphs yield quantities that fall as exp[— E(\T\ + |T"|)]. The 
matrix element is given by the limit as T, T' — > oo of the ratio of the upper graph to the 
lower graph, with the same choice of sources and sinks in both graphs, times the coefficient 
of the exponential fall off for the point-point quarkonium propagator. We used noisy-point 
and noisy-gaussian sources and generated retarded and advanced quark propagators from 
each time slice. We chose the Coulomb gauge for the field configurations. This choice made 
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implementation of extended sources simpler and allowed us to replace covariant derivatives 
with normal derivatives, with errors of relative order v 2 . We calculated heavy-quark propa- 
gators (j?(x, t) on the lattice, using the nonrelativistic formulation of Lepage et al. ||, with 
an evolution equation that is valid to the lowest non-trivial order in v 2 : 

G{x,t+ 1) = (1 - H /2n) n Ul t {l - H Q /2n) n G(x,t) 

+ 5x,(A+i,o, (4) 

with G(x, t) = for t < 0, and H = —A^/2M — h . A^ is the gauge-covariant discrete 
Laplacian, M the bare heavy-quark mass, and h = 3(1 — u )/M , with u = (\U p i a qette)^ ■ 
We chose n = 2. Since Ti/Mq is suppressed by 0(v 2 ) relative to Qi, it is of the same order 
as terms that we have neglected in the computation of Q\. The main justification for its 
calculation is that there are decays such as 3 Si — > Light Hadrons, 3 Si — > 7 + Light Hadrons, 
and 3 Si — > 37, for which the coefficient of Fi/Mq is approximately —5 times that of Qi [§[0. 
In these cases, the contributions of terms proportional to T\ could be important. 

In our lattice calculations we used 149 quenched gauge-field configurations on a 16 3 x 
32 lattice with 6/g 2 = 6.0 for bottomonium and 158 configurations on an 16 3 x 32 lattice 
with 6/g 2 = 5.7 for bottomonium and charmonium. For bottomonium we took M = 1.5 
at 6/g 2 = 6.0 and M = 2.7 at 6/g 2 = 5.7. For charmonium at 6/g 2 = 5.7 we took 
Mo = 0.69. These values correspond to those used by the NRQCD collaboration ||. (Note 
that their mass definitions are u times ours.) The values of u that we used are 0.87778701 
at 6/g 2 = 6.0 and 0.8608261760 at 6/g 2 = 5.7. Except where we explicitly state otherwise, all 
quantities in this paper are in lattice units. To convert to physical units, we use inverse lattice 
spacings a~ l = 2.4 GeV for bottomonium at 6/g 2 = 6.0, a~ l = 1.37 GeV for bottomonium 
at 6/g 2 = 5.7, and a" 1 = 1.23 GeV for charmonium at 6/g 2 = 5.7. These are the values 
obtained by the NRQCD collaboration. Our error estimates do not include the errors in 
these quantities. 

NRQCD predicts that 

3 



^ 1 /l( 1 5o|^ t x|0)| 2 = (l + ^ 4 )) (5) 
Wi/|( 1 P 1 |^YDx|0)r = (l + ^ 4 )), (6) 

where the vacuum-saturation approximation amounts, in this case, to ignoring the 0(v A ) 
term. For bottomonium at Q/g 2 = 6.0 we measured the v 4 term for Q\ to be 1.3(1) x 1CT 3 . For 
charmonium at Q/g 2 = 5.7 this term is approximately 1%. For TCx, these (9(v 4 ) terms, while 
larger than those for Qx, are still quite small. Thus, the vacuum saturation approximation 
is even better than one would expect. We will therefore use the vacuum-saturation values 
for Qi, Hi, and T\ in the discussions to follow. The lattice quantities Q\l and TCn are then 
given by the coefficients of the exponentials in the S- and P-wave quarkonium propagators, 
respectively. 





charmonium 


bottomonium 


Q/g 2 


5.7 


5.7 


6.0 


Gil 


0.1317(2)(12) 


0.9156(9)(65) 


0.1489(5)(12) 


T 1L {non)/G\L 


1.2543(7) 


2.7456(8) 


1.3135(8) 


T 1L {cov) 1 VlL 


0.5950(5) 


2.1547(7) 


0.8522(5) 


J r 1L (non 2 ) i 'Gil 


0.7534(4) 


1.2205(2) 


0.7775(5) 


Fil{cov 2 ) / Gil 


0.5201(3) 


1.1111(2) 


0.6659(3) 


T~tlL 


0.0208(2)(20) 




0.0145(6)(20) 


'Hsl/'Hil 


0.034(2)(8) 




0.0152(3)(20) 



TABLE I. Lattice decay matrix elements expressed in lattice units (a = 1). Note that P-wave 
bottomonium matrix elements have yet to be calculated at Q/g 2 = 5.7. 



A summary of our results for the lattice matrix elements defined in (|3|) is presented 
in Table |. When a second error has been included, it is an estimate of the systematic errors 
associated with the parametrization of the fitting functions and with the contamination from 
higher states for propagators in which the separation between source and sink is too small. 

To the order in v 2 in which we are working, our lattice matrix elements are related 
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to their continuum counterparts by 



Gil 



(7a) 



Til 



(1 +7)^1 + ^1, 



(7b) 



and 



HiL 



(1 + l)Hi + kH s , 



(8a) 



^8L 



(1 + rj)Hs + (Hi 



(8b) 



where the subscript L indicates the lattice quantity. The coefficients e, 7, 0, t, 77 and £ 
are of order a s ; k is of order a 2 . We have calculated these coefficients through order a s 
(one loop) in tadpole-improved perturbation theory ||. Our values for these coefficients, for 
MS regularization of the continuum matrix elements, are given in Table ||. The accuracy of 
the coefficients of a s in this table is estimated to be better than 1%. In computing (, we have 
taken the factorization scale to be 1.3 GeV for charmonium and 4.3 GeV for bottomonium. 
These values correspond, approximately, to the MS heavy-quark masses. Note that </> and k, 
in physical units, have dimensions of (mass) 2 , and ( has dimensions of 1/mass 2 , whereas the 
other coefficients are dimensionless. If we render 0, k, and ( dimensionless by dividing T\ 
and (f) by Mq, by dividing Hi and k by Mq, and by multiplying ( by Mq, respectively, then 
none of the coefficients of a s is exceptionally large. Hence, the use of low-order perturbation 
theory appears to be reasonable. 

For those coefficients that arise from a positive integrand, the method of Lepage 
and Mackenzie |J yields an optimal scale for a s that is close to 1/a. Thus, we choose 
a s = a v (l/a) = 0.3552 at 6/g 2 = 5.7 and 0.2467 at 6/g 2 = 6.0. 

Substituting the numerical values from Tables | and |I| into (0) and @, we obtain the 
results shown in Table [TIT]. In Table |TTT], the first and second errors in the lattice results are 
from the statistical and systematic errors in Table |. The third error is an estimate of the 
systematic error that arises from the neglect of terms of higher order in a s in the coefficients 
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charmonium 


bottomonium 


6/g 2 


5.7 


5.7 


6.0 


e 


-0.7326 a s 


0.2983 a s 


-0.4877 a s 


~f(non) 


-0.02578 a s 


-1.248 a s 


-0.9117 a s 


~f(cov) 


-2.860 a s 


-2.192 a s 


-2.560 a s 


■y(non 2 ) 


-0.2774 a s 


-1.096 a s 


-0.9236 a s 


4>(non) 


1.486 a s 


10.90 q s 


4.418 a s 


4>(cov) 


0.3928 a s 


9.808 a s 


3.325 q s 


<p{nori2) 


1.004 a s 


6.096 q s 


2.863 q s 


l 


-0.7603 a s 


-1.852 a s 


-1.191 a s 


V 


0.09157 a s 


-0.03728 a s 


0.06096 a s 


c 


-0.1785 a s 


-0.006011 a s 


-0.01862 a s 



TABLE II. Coefficients relating lattice and continuum matrix elements. The arguments of 7 
and (ft correspond to different lattice representations of T\. cov is a tadpole-improved Q naive 
discretization of the gauge covariant continuum operator; non is the simple, gauge- noncovariant, 
finite-difference operator in coulomb gauge; the subscript 2 indicates a difference operator with 
spacings of 2 lattice units. 



of Table ||. It is obtained by taking the uncertainty in the coefficients to be either ot 2 s times 
the zeroth order term (if any) or a s times the magnitude of the first order term, whichever 
is the larger. In the case of J-'i/Qi, the uncertainty is large, so we have presented our results 
as a range of values. 

For purposes of comparison, we have also shown in Table |II| the experimental (phe- 
nomenological) results, where available, for the matrix elements that we have computed. 
The phenomenological results for Q\ were extracted from the measured decay rates for 
J/i[) —> e + e~, T] c — > 77 and T — > e + e~ |7j], using the expressions in Ref. Values for 
T\jQ\ for Jjijj are those of Ko, Lee and Song ||. The results for Hi and Tig/Hi for Xc are 
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lattice 


experiment 




lattice units 


physical units 




charmonium 6/g 2 = 5.7 






Qi 


0.1780(3) (16) t^t 


0.3312(6) (30) till Gey3 


0.36(3) GeV 3 




0.05 — 0.54 


0.07 — 0.82 GeV 2 


0.057 GeV 2 




0.0285(2) (27) 


0.0802(6) (77)t\il Gey5 


0.077(19)(28) 
GeV 5 


Hs/Hi 


0.086(l)(6)ii 


0.057(1)(4) + 5T GeV" 2 


0.095(31)(34) 
GeV" 2 


bottomonium 6/g 2 = 5.7 






Gi 


0.8279(8) (59)1^6 


2.129(2)(15)1 2 ^ GeV 3 


3.55(8) GeV 3 




-3.7 — 0.2 


-6.9 — 0.4 GeV 2 





bottomonium 6/g 2 = 6.0 






Qi 


0.1692(6)(14)±gg 


2.340(8) (19)t\ll GeV 3 


3.55(8) GeV 3 


Fi/Gx 


-0.34 — 0.28 


-2.0 — 1.6 GeV 2 




Hi 


0.0205(9) (28) tig 


1.63(7)(23)±i§ GeV 5 




Hs/Hi 


0.0151(2) (14) 


0.00262(3)(24)l^GeV- 2 





TABLE III. Continuum MS decay matrix elements from our lattice calculations, compared 
with those extracted from experimental decay rates, where available. 



from ||; the first error is experimental, the second theoretical. For P-wave bottomonium, 
there is as yet no published data on decays into light hadrons, photons, and/or leptons. 
The extraction of phenomenological matrix elements from the experimental data requires 
values for the heavy-quark masses. Our choices correspond to pole masses of 5.0 GeV for the 
b quark, the result obtained by the NRQCD collaboration [Q, and 1.5 GeV for the c quark 
IJTIJ. We also require values for a s and a in order to evaluate the partonic decay rates. For 
these we used a s {M c ) = 0.243, a s {M b ) = 0.179, a(M c ) = 1/133.3 and a(M b ) = 1/132. 
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In the above analysis, we have not taken into account the errors due to the omission 
of terms of higher order in v 2 . These could be as large as 10% for bottomonium and 30% 
for charmonium. For Gil, the NRQCD collaboration has published results that are accurate 
to the next-to- leading order in v 2 . Since these higher-order results distinguish the singlet 
and triplet states, we compare our results with a weighted average of their results. For 
charmonium at Q/g 2 = 5.7 they obtain Gil = 0.133(4), and for bottomonium at Q/g 2 = 6.0 
they obtain Gil = 0.144(4), in good agreement with the results of Table |. For these matrix 
elements, as with masses, most of the effect of contributions of higher order in v 2 is to split 
the results for the singlet and triplet states without shifting the weighted average. 

There are some additional sources of error that we have not included in Table fH\ . 
One of these is the uncertainty in the physical value of a -1 . Using the results of the NRQCD 
collaboration, we find that the uncertainties in the values of the matrix elements from this 
source are 7% for G\ and 13% for 7i x in charmonium and 13% for Gi and 23% for TCi in 
bottomonium. As we have already mentioned, extraction of the phenomenological matrix 
elements requires knowledge of the heavy-quark mass. The NRQCD collaboration quotes an 
error of 4% for the 6-quark mass, which introduces an 8% error in Gi and a 16% error in TCi- 
For the c-quark mass we have no good error estimate. Further sources of uncertainty are 
the QCD radiative corrections to parton-level decay rates. Estimates of these uncertainties 
have been included in the phenomenological values of 7i\ and 7i% for charmonium that are 
reported in Ref. @. Finally, there are the errors that arise from using quenched (rather 
than full) QCD, for which we have no estimates. 

Let us now discuss our results. For charmonium Gi, Hi, and Hg/Hi are in agreement 
with experiment, although both the lattice and experimental results have sizeable errors. We 
note that we would not have found this agreement had we failed to include the perturbative 
corrections that relate the lattice matrix elements to the continuum ones. The quantity 
Fi/Gi is poorly determined, for both charmonium and bottomonium, owing to the mixing 
of T\ with Qi in (0). Because T\j{M 2 G\) is of order t) 2 C 1, a coefficient 0/M 2 of order 
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a s yields a large mixing, and any uncertainties in 0/M 2 are amplified in T\j [M 2 Q\). We 
do learn, though, that J r i/(M 2 Q 1 ) is no larger than 0(v 2 ), in agreement with the NRQCD 
scaling rules 0,[]J. For charmonium, T\jQ\ is probably positive, while for bottomonium a 
negative value is preferred. In the case of bottomonium, the lattice result for Q\ is 35 - 
40% below the experimental value, although there is good agreement between the 6/g 2 = 5.7 
and 6/g 2 = 6.0 predictions. At least part of this discrepancy, which was first noted by the 
NRQCD collaboration, is due to the quenched approximation Our results for the P- 

wave matrix elements for bottomonium can be translated into predictions for bottomonium 
decay rates ||. In the P-wave case, these should lead to significant new tests of the theory 
as the relevant experimental data become available. 

As is clear from Table [TTJ, the largest uncertainties in the lattice matrix elements 
(aside from those due to quenching) come from neglecting higher-order corrections to the 
coefficients of Table [TT]. This suggests that a useful strategy might be to use lattice methods 



HI to compute the coefficients beyond leading order. In addition, one might consider the 



use of alternatives to the MS regularization of the continuum matrix elements so as to 
avoid renormalon ambiguities in the matrix elements and short- distance coefficients |TT 



Note that, to the extent that we can replace Tin and 7i\ by their vacuum saturation 
approximations, Tii and 7i% are free of renormalon ambiguities. Of course, these ambiguities 
cancel in physical quantities if one works consistently to a given order in a s in both the 
lattice-to-continuum coefficients and the short-distance coefficients. 

It is an interesting fact that, for decays of P-wave states in both charmonium and 
bottomonium, the ratio of the octet matrix element to the singlet matrix element is in 
reasonable agreement with a crude phenomenology. This phenomenology is based on solving 
the one- loop evolution equation for 7i 8 (Ref. |T|) and assuming that 7Yg vanishes below a scale 
Mqv. Since the one-loop evolution contribution to the decay matrix element 7ig is the same 
as that for the corresponding production matrix element 7i' s , this simple phenomenology 
suggests that TC' 8 is approximately equal to 7^8- For charmonium production, TC' 8 has been 
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extracted from CDF data [|1^] and can also be deduced from recent CLEO data |14| . Dividing 
the CDF and CLEO values of 7ig by the phenomenological value of Tlx given in Table [TTT], we 
obtain 0.042(19)GeV -2 and 0.046(28)GeV -2 respectively, which agree, within errors, with 
our lattice result for TCsfHi- 
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